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Using particle-in-cell simulations, we examine hot electron generation from electron plasma waves 
excited by stimulated Raman scattering and rescattering in the kinetic regime where the wavenumber 
times the Debye length {kXn) is ^ 0.3 for backscatter. We find that for laser and plasma conditions 
of possible relevance to experiments at the National Ignition Facility (NIF), anomalously energetic 
electrons can be produced through the interaction of a discrete spectrum of plasma waves generated 
from SRS (back and forward scatter), rescatter, and the Langmuir decay of the rescatter-generated 
plasma waves. Electrons are bootstrapped in energy as they propagate into plasma waves with 
progressively higher phase velocities. 
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Stimulated Raman scattering (SRS), the decay of a 
light wave into a forward propagating electron plasma 
wave (EPW) and a forward propagating (SRFS) or back- 
ward propagating (SRBS) light wave, involves fundamen- 
tal nonlinear wave- wave and wave-particle interactions. 
SRS continues to be studied extensively because the loss 
of incoming energy due to backscatter and the potential 
fuel preheat due to hot electrons generated by the EPW 
are threats to Inertial Fusion Energy (IFE) devices such 
as the National Ignition Facility (NIF). Recent experi- 
ments at NIF have shown electron heating up to energies 
above 100 keV [1 . A low-temperature (Tg = 10-20 keV) 
part of the heated electron distribution can be attributed 
to SRBS, but the high-temperature part is currently un- 
explained. There is speculation that these electrons are 
generated near the quarter-critical density by instabili- 
ties such as two-plasmon decay or SRFS [2^. 

In this article, we present a novel mechanism for gen- 
erating 100 keV electrons through SRS rescatter, specif- 
ically through SRBS of SRBS, SRBS of SRFS, and the 
Langmuir decay instability (LDI) of the rescatter plasma 
waves, where LDI is the decay of an EPW into a counter- 
propagating EPW and an ion acoustic wave. We further 
show how electrons can get progressively heated as they 
travel between waves of increasing phase velocities. This 
mechanism allows rescatter and SRFS to heat electrons 
that have been initially heated by SRBS, even though 
the phase velocity of SRFS is too high for it to trap and 
heat electrons on its own. 

Particle-in-cell (PIC) simulations have been used to 
study rescatter and multi-stage electron heating from 
SRS, albeit only electron heating between SRBS and 
SRFS. Hinkel et al. [3 showed rescattering for NIF- 
relevant parameters but not the resulting hot electrons. 
Other authors [IHl] have shown electron heating by 
SRFS, in some cases explicitly due to SRFS accelerating 
electrons initially heated by SRBS, but these simulations 
have been for more intensely-driven and/or hotter plas- 
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FIG. 1: Electron phasespace showing trapped particles dur- 
ing SRBS (left) with the corresponding flattening of the distri- 
bution function (right, dashed) and initial distribution (right, 
solid). 



mas outside the current range for NIF where electron 
temperatures Tg ~ 2-6 keV and laser intensities / and 
wavelengths Aq are such that /Aq ~ 10^^ W /im^/cm^ 
(in laser hot spots). 

The electron energies that result from trapped parti- 
cle interactions in an EPW depend on its phase velocity, 
and potential amplitude, In the wave frame, the 
energy of a plasma electron, (7' — l)mc^ — e<I>', is con- 
served, where ^' = and 7 = . =. A trapped 

electron with the highest energy will be at rest at the top 
of the potential hill in the wave frame. At the bottom 
of the potential such an electron has a 7^^^ = 1 + 

and = ±x/(l + ^^)^ - 1, where A$ is the po- 

tential difference between the top and bottom of the 
wave. 7max, Vmax, and Vmax cau then be obtained by 
transforming these expressions back to the lab frame. In 
the nonrelativistic limit, the maximum velocity and en- 
ergy of a trapped electron are Vmax = '^0 + Y^2e(A<I>)/m 
and Eraax — ^^^^max' ^ siuusoidal plasma wave, 

= 2^rnax and Vmax = ^0 + ^tr, with the trapping 
width Vtr = 2yfe^Jm. Figure [l] illustrates trapped elec- 
trons oscillating between Vmax = '^0 ± Vtr in a kinetic 
SRS simulation (the first simulation detailed later). 
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EmaxforTe = 2.5keV 

^ J — nonrelativistic 
> 200- ^relativistic 
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FIG. 2: Estimate of bound on maximum electron veloc- 
ity Vmax (top-left) and kinetic energy Emax (assuming Te 
— 2.5 keV, top-right). EPW phase velocities for one set of 
parameters (bottom) illustrate the phase- velocity-ordering of 
modes, with rescatter and LDI of rescatter intermediate be- 
tween SRFS and SRBS. 



One can also estimate a bound on maximum electron 
energy by estimating that the plasma wave is at the warm 
wavebreaking limit |8]. The wavebreaking derivation of 
[9] shows that extrema in <l>, 
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occur for roots {E = 0) of: 
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where a waterbag distribution is assumed, E = 
eE/mujpV(f) is the normalized electric field, v = v/v^j^^ 
/3 = 3{vth/V(f))'^ ^ and cq is an arbitrary constant. Assum- 
ing Epeak is the wavebreaking amplitude, two of the roots 
in V are approximately evenly spaced about and rep- 
resent fluid velocities at the extrema in ^. We can then 
use these two roots to calculate A$ between the extrema 
and substitute it into the above expressions for maximum 
energy. Figure [2] top- left shows Vmax/'^cf) as a function of 
/3, from which it is seen that for /3 > 0.1 {kX^ ^ 0.18) 
the difference between Vmax and is no bigger than 'u^, 
i.e., Vmax < ^v^. In addition, A<l> = 2^rnax for /3 > 0.1 
since the plasma wave is nearly sinusoidal. Figure [2] top- 
right shows £max assumiug that Tg = 2.5 keV, where the 
dotted line is from £max = \^^max solid line 

includes relativistic corrections. In simulations, we find 
that Epeak i^peak) fo^ the SRBS wave is typically ^2/3 
of the wavebreaking estimate, so these curves should be 
viewed as a limit. 

The appropriate and (kinetic energy for a parti- 
cle at ^0) for the various plasma waves are shown in Fig- 
ure [2]-bottom. The plots in combination show that SRBS 
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FIG. 3: Top) Temporal evolution of the SRS EPW wavenum- 
bers. Bottom) Electron energy spectrum. 



does not longitudinally accelerate electrons to 100 keV ki- 
netic energies. Trapped particles with additional trans- 
verse velocity components may reach such energies, as 
may have occurred in L. Yin et al. pTl, but we leave this 
for future work. For rescatter, on the other hand, with 
higher (and lower /3), 100 keV is well within range. 
Rescatter EPWs, as well as their LDI decay EPWs, are 
potential producers of 100 keV electrons. 

In this article, we simulate the scattering processes in 
one and two dimensions (ID and 2D) using the electro- 
magnetic PIC code OSIRIS [10]. The electrons have a 
temperature Tg = 2.5 keV (3 keV) in ID (2D) and density 
between Ue = 0.09 — 0.1 Oner {kXo ^ 0.33 for backward 
SRS); ions are either fixed or mobile with ZT^/Ti = 2 
and Mi/rrie = 1836. The laser has normalized peak 
electric field eE/mcuJo = 0.0164 which corresponds to 
/o = 3 X 10^^ W/cm^ for Aq = 0.351/im. In 2D the laser 
is focused from /q = 3 x 10^^ W/cm^ at the simulation 
edge to /q = 5 X 10^^ W/cm^ at focus with a focal spot 
size of 2.6 /im (8Ao). The laser propagates along x and 
is polarized in z. The simulations have absorbing bound- 
aries for the fields and thermal-bath boundaries for the 
particles, with 16384 (16384 x 512) cells and 512 (256) 
particles per cell in ID (2D) to simulate plasma of size 
180 jam (200 x 15 jam^). The length corresponds to an 
f/S speckle of length S fXo = 180/im. 

First we consider a ID run with immobile ions and 
homogenous density Ue/ricr = 0.10. The temporal evo- 
lution of the SRS plasma wave wavenumbers can be seen 
in Figure [3| along with the corresponding temporal evo- 
lution of the electron distribution function. SRBS grows 
first, as it has the largest growth rate. The growth of 
SRFS follows, and after SRFS has grown to sufficient 
amplitude for its daughter light wave to be above the 
SRS threshold, the SRFS saturates via SRBS of SRFS. 

The hot electron tails in the distribution follow a dif- 
ferent progression. The daughter plasma wave phase 
momenta {pe/rrieC = jv^p) increase from SRBS (0.26) 
to SRBS of SRFS (0.56) to SRFS (2.4), with electrons 
at those speeds having kinetic energies of 17, 75, and 
820 keV respectively. Electron trapping by SRBS starts 
at ujot ~ 10000; this process does not accelerate elec- 
trons above 70 keV. SRFS grows to a mode amplitude 
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larger than SRBS by ujot = 20000, but normalized to its 
wavebreaking value it is smaller so it does not trap par- 
ticles and has no immediate effect on the hot electron 
tail. Electrons begin to be accelerated to energies above 
70 keV by the rescatter that develops at uJot ~ 20000, 
and by uJot ~ 33000 electrons have been accelerated to 
sufficient energies by the rescatter that SRFS can inter- 
act with a significant number of electrons, trapping them 
and accelerating them beyond 250 keV all the way up to 
1 MeV. With heating by both rescatter and SRFS, ap- 
proximately 0.1% of the electrons get heated above 100 
keV. 

Though not shown, the maximum EPW amplitude in 
the region of SRFS activity is eE/mV(f)UJp ^ e^/mv'^ ^ 
0.37 at UJot ~ 33000. Using = 0.93c, the energy 
an electron needs in order to be trapped is 140 keV 
{v = 0.62c). Since SRBS only generates electrons with 
energies less than approximately 60 keV, this illustrates 
why SRFS requires the intermediate step of rescatter. 
This is consistent with Figure |2] where the electron en- 
ergy sweeps to higher energies once rescatter heats them 
to 140 keV. 

Rescatter can be limited by several effects. We have 
performed several simulations to study its dependence 
on density. Even for strong SRS, rescatter does not 
grow above the scattered light's quarter-critical density, 
which for scattered light of frequency uj ^ ujq — ujp is 
n/ucr ^ 0.11. On the other hand, for lower densities 
such as n/ricr < 0.09, the growth rates of all SRS pro- 
cesses decrease, likewise making rescatter less likely. In 
addition, both instabilities are sensitive to density gra- 
dients and other scattering instabilities such as LDI or 
Brillouin scattering of the scattered light. 

While we do not see Brillouin scattering for our pa- 
rameters, we do see saturation by LDI and quenching of 
SRFS by density gradients. Here we turn to a ID mo- 
bile ion simulation with a linear density gradient from 
Tie/ Tier = 0.09 to 0.10 ovcr the length of the domain (180 
/im). In this simulation, SRFS does not grow. Never- 
theless, SRBS and SRBS of SRBS light stih occur, as 
does LDI of the rescatter plasma wave. The spectrum of 
plasma modes can be seen in Figure |4[top. 

Figure [4]-bottom shows that all of the plasma waves 
heat electrons, with rescatter (here SRBS of SRBS) again 
accelerating electrons up to energies of 100-200 keV. Fur- 
thermore, the EPW from LDI decay of rescatter also 
heats electrons, and as it travels in the opposite direction 
as the rescatter EPW, the combined instabilities gener- 
ate energetic electrons in both directions. The LDI de- 
cay product has a slightly lower wavenumber compared 
to the decaying EPW, interacting with the electron dis- 
tribution at slightly higher phase velocities. The hot tail 
due to LDI therefore extends to higher energies than the 
tail due to rescatter. Furthermore, even though the LDI 
EPW interacts with the electron distribution farther out 
in its tail, it heats more total electrons than the rescatter 
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FIG. 4: Top) frequency vs wavenumber of El-field; bottom) 
electron distribution versus electron kinetic energy for /o = 3 
(solid) and 4 (dashed) xlO^^ W/cm^, where +/- represent 
forward/backward traveling electrons. 



EPW since it also interacts with the previously formed 
hot tail from the original SRBS. 

If we apply our earlier theoretical estimate of electron 
^max to the rescatter and LDI plasma waves, we can 
test this limit against the electron spectrum shown in 
the Figure. The phase velocities can be ascertained from 
the spectrum, where the EPWs have ~ for 
rescatter of SRBS and ^Avth for LDI of rescatter. The 
electron kinetic energies corresponding to the theoreti- 
cal limit (Fig. [2| of Vmax are ^ 200 keV for SRBS and 
300 keV for LDI. These indeed bound the upper edges 
of the flat tails in the plotted distribution in Figure [4j 
even for an exactly similar case with higher laser inten- 
sity. The low-velocity end of the hot tails correspond 
approximately to ^ mc^(70 — 1) (?^ 44 and 61 keV). 

Finally, we present results from a 2D simulation of 
a single speckle. Figure [5}top-left shows the temporal 
evolution of wavenumbers for the Ez (transverse) field 
along the central axis. The bursty mode at kc/uj^ ^ 0.5 
corresponds to SRBS, while the steadily growing mode 
at kc/ujQ ~ 0.6 that peaks at uo^t ~ 17000 corre- 
sponds to SRFS (the anti-Stokes mode is also present at 
kc/ujQ ~ 1.3). Rescatter of both light waves is present, 
with SRBS of SRBS at kc/ujo ~ 0.20 and SRBS of SRFS 
at kc/ujQ ^ 0.15 starting at uJot ~ 17000. Correspond- 
ing rescatter plasma wave modes are seen in the Ex field 
(not shown), as well as broadband signals from modes 
that grow after ujot ~ 17000 due to LDI. 

The electron distribution shown in Figure [5}top-right 
fiattens slightly at uJot = 10000 due to the first burst of 
SRBS, followed by much more energetic tails at uJot = 
18000 due to rescatter and LDI of rescatter. Though not 
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FIG. 5: Top) Temporal evolution of transverse E- field 
wavenumbers along the simulation center (top-left) and lon- 
gitudinal electron distribution spatially averaged over x and 
y (top-right). Bottom) Charge density of electrons at Vcf) of 
BSRS (grey), SRBS of SRFS (blue), and SRFS (orange). 



Electron Distribution Functions: 
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FIG. 6: Energy distributions. The 2D case includes trans- 
verse momentum, with dz referring to the sign of px- 



shown here, the electron phasespace reveals that the pos- 
itive momentum tail is caused by SRBS of SRFS and the 
negative momentum tail by LDI of SRBS of SRFS. The 
importance of trapped electron bootstrapping between 
SRBS and rescatter can be seen in the bottom plot of 
Figure [5) where we plot the charge density amplitude in 
electron phasespace as a function of space and time for 
three different ranges of electron momenta. The phases- 
pace bins Pe/rrieC = (0.32, 0.56), (0.56, 0.72), and (0.76, 
2.50) cover V(fy of the plasma waves due to SRBS, SRBS of 
SRFS, and SRFS respectively. With SRBS growing be- 
hind the laser focus (focus at xloq/c = 1790) and SRFS 
growing in front of the laser focus, the electrons trapped 
and accelerated by SRBS first have to cross the simu- 
lation length before interacting with the region where 
SRFS (and rescatter of SRFS) has grown. After they do 
so (as shown in grey, the rescatter can interact with these 
electrons and accelerate them further. The blue color 
represents electrons heated by rescatter once the elec- 
trons heated by SRBS enter the region of rescatter, while 
the orange color shows further acceleration by SRFS. LDI 
limits SRFS for uot > 20000, and thereby also rescatter 
of SRFS. At uot ^ 18000, those hot electrons with kinetic 
energies above 100 keV have a forward- going kinetic en- 
ergy flux of approximately 3% of the total incident laser 
poynting flux, while subsequent fluxes at ujot ~ 36000 
and 43000 are both ^ 0.2%. 

The electron distributions in energy for the ID immo- 
bile ion case and the 2D case are shown in Figure [6j where 
one can see that electrons are not accelerated to 100+ 
keV energies until rescatter has grown {ooot > 17000). 
Fitted lines for temperatures show that the range of elec- 
tron energies, not the slope of the distribution, identifies 
which plasma wave (instability) is responsible for those 



hot electrons. 

The range of energies shown in this article is consistent 
with reported results of hot electron measurements from 
NIF and shows that SRS rescatter should be considered 
as a source of 100 keV electrons. While the results here 
are limited to single speckle physics with laser intensities 
at the higher end of hot spot intensities, one might rea- 
sonably assume that scattered light will be amplified to 
levels seen here as it travels through multiple speckles. 
Yin et al. [11 have shown 100 keV electrons in multi- 
speckle SRS, possibly resulting from electrons interacting 
with multiple SRBS plasma waves; however their density 
parameters were outside the allowable range for rescat- 
ter. Studying rescatter in multi-speckle simulations is an 
area for future work. Furthermore, Yan et al. [12 have 
shown multi-stage electron acceleration in two-plasmon- 
decay simulations where there occurred a broad range of 
plasma waves (in direction and phase velocity). We an- 
ticipate that the bootstrap acceleration of hot electrons 
traveling between speckles and separate regions of SRS 
is worth further study, as well as the intriguing idea of 
whether scattered light can undergo not just SRS rescat- 
ter but also two-plasmon decay as it travels from regions 
above its quarter-critical density to lower densities. 
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